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ABSTRACT 


Plane  strain  fracture  toughness  and  tensile  properties  were  determined  at  room  tempera¬ 
ture,  utilizing  compliance  and  pop-in  methods,  for  four  high-strength  steel  alloys:  PH  15-7 
Mo,  17-7  PH,  AM  350  and  Vasco  Jet  1000  (H-11).  Fracture  toughness  values  varied  over  a 
fairly  wide  range,  with  AM  350  having  the  highest  at  approj^ately  60.7  KSI-v/in.  and  Vasco 
Jet  1000  (H-11)  having  the  lowest  at  approximately  25  KSIvin.  A  computer  program  used  to 
reduce  fracture  toughness  data  was  able  to  calculate  critical  crack  length  as  well  as  fracture 
toughness  when  given  either  suitable  compliance  gage  data  or  the  measured  test  data.  An 
acoustical  pickup,  used  as  an  additional  test  monitor,  is  described.  Analytical  basis  for  the 
compliance  method  is  presented. 
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INTRODUCTION 

As  the  quest  for  better  performance  aerospace  systems  continues,  material  strength  re¬ 
quirements  become  higher  and  higher.  The  use  of  impact  energy  absorption  test  data  has 
become  inadequate  for  describing  the  toughness  characteristics  of  materials.  In  the  continuing 
search  for  high-strength  material,  it  becomes  more  and  more  important  to  know  (with  ac¬ 
curacy)  thebehaviorof  a  material  (as  a  formed  aerospace  component)  under  service  conditions. 
Therefore,  a  proposed  material  must  be  evaluated  for  its  resistance  to  brittle  fracture  in  the 
thicknesses,  strength  levels,  and  metallurgical  conditions  envisioned.  If  such  an  evaluation  is 
not  made,  the  behavior  of  the  material  is  not  completely  characterized  and  either  overdesign 
or  underdesign  can  result.  Both  cases  are  undesirable  —  overdesign  will  result  in  a  heavy, 
inefficient  structure,  while  underdesign  can  result  in  catastrophic  failure. 

The  plane  strain  fracture  toughness  parameter,  gives  to  the  designer  some  measure 

of  service  performance  when  a  small  flaw  is  present  in  the  material.  By  knowing  the  size  of 
the  largest  flaw  which  can  be  tolerated  in  a  structure,  the  designer  can  determine  the  stress 
level  which  will  cause  failure,  and  then  design  with  an  appropriate  safety  factor.  Alternatively, 
he  can  determine  the  smallest  flaw  detectable  with  available  inspection  equipment,  determine 
the  stress  level  which  will  cause  failure  with  this  flaw  present,  and  design  below  that  stress. 

Four  high-strength  steel  alloys  in  sheet  form  were  tested  in  this  program.  This  report 
shows  that  the  plane  strain  fracture  toughness  of  high-strength  alloys  is  very  important  even 
for  relatively  thin  sheet. 


MATERIALS 

Table  1  shows  a  list  of  the  four  sheet  materials,  their  thicknesses,  and  heat-treatment  con¬ 
ditions  used  in  this  program.  Standard  heat-treatment  designations  were:  PH  15-7  Mo,  RH  950; 
17-7  PH,  RH  1050;  AM  350,  SCT  950;  and  Vasco  Jet  (VJ)  1000  (modified  AISI  type  H-11), 
triple  tempered  at  1000°F.  The  alloy  suppliers,  heat  numbers,  and  the  chemical  composition 
of  the  materials  used  in  this  program  are  also  listed.  Photomicrographs  of  the  four  sheet 
materials  in  the  heat-treated  condition  are  shown  in  Figure  1. 


SPECIMENS 

The  fracture  toughness  tests  conducted  in  this  program  used  a  3  in.  by  12  in.  center  notch 
specimen.  A  drawing  of  this  specimen  is  shown  in  Figure  2. 

The  specimen  preparation  procedure  was  as  follows:  shear  to  size,  drill  1  in.  dia  pinholes, 
heat-treat  in  fixtures,  surface  grind  to  remove  scale,  “Elox”  center  notch,  drill  shim  holes, 
and  tension-tension  fatigue  crack  the  starter  notch  to  approximately  1  in.  (1/3  W).  The  center 
notch  was  placed  in  the  specimen  after  heat  treatment  to  avoid  any  he  at- treatment  distortion 
problems.  The  specimens  were  heat-treated  in  an  atmosphere  of  dissociated  ammonia. 

Recently,  the  ASTM  committee  on  Fracture  Testing  of  High-Strength  Metallic  Materials 
recommended  that  the  specimen  width-to-thickness  ratio  should  be  between  16  and  45.  How¬ 
ever,  the  specimens  for  this  program  were  machined  prior  to  this  recommendation  and  the 
width-to-thickness  ratios  ranged  from  50  for  VJ  1000  to  85  for  17-7  PH.  Because  of  the  rel¬ 
atively  thin  sheet  material  used,  considerable  difficulty  was  experienced  with  buckling  around 
the  pinholes.  To  prevent  this  buckling,  hardened  steel  carrier  shims  were  used.  These  shims 
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were  fastened  to  the  specimen  by  pins  located  in  the  three  peripheral  holes  shown  in  Figure  2, 
These  carrier  shims,  while  not  fastened  to  the  grips,  created  more  load  bearing  surface  and 
distributed  the  stress  around  the  pinholes.  This  method  of  shimming,  as  described  in  the  first 
ASTM  Fracture  Toughness  of  High-Strength  Steel  Committee  Report  (Reference  1)  offers  a 
fast  and  convenient  means  to  eliminate  pinhole  buckling  problems. 

Standard  2-in.  gage  length  smooth  tensile  specimens  were  machined  and  heat-treated  si¬ 
multaneously  with  the  fracture  toughness  specimens. 


TEST  EQUIPMENT 

All  tensile  and  fracture-toughness  testing  was  done  on  a  50,000  lb  Baldwin  Universal  tensile 
machine.  The  fracture  specimens  were  fatigued-cracked  in  tension- tension  using  a  6-ton 
Schenck  fatigue  machine  shown  in  Figure  3.  Fatiguing  at  a  maximum  stress  of  1/5  the  yield 
stress  at  2000  cycles/min  produced  fatigue  cracks  of  the  desired  length  in  approximately 
15  minutes.  The  brittle  behavior  of  VJ  1000  necessitated  a  lowering  of  the  maximum  stress 
after  crack  initiation  to  prevent  the  fatigue  crack  from  becoming  unstable  as  it  approached 
the  desired  length  of  1  in.  Fatigue  crack  lengths  were  measured  after  testing  with  a  toolmaker’s 
microscope. 

The  compliance  gage  chosen  was  originally  developed  by  Richard  W.  Boyle  at  the  Naval 
Research  Laboratory  (Reference  2).  This  compliance  gage  was  modified  to  use  a  Model  PS-6M 
Baldwin  Microformer  which  provides  a  total  magnification  of  approximately  550X.  The  com¬ 
pliance  record  was  made  on  a  Baldwin  autographic  recorder. 

During  the  program  it  was  found  that  the  magnification  produced  by  this  combination  of 
microformer  and  gage  was  insufficient  to  reliably  record  the  pop-in  event.  It  was  then  decided 
to  obtain  an  acoustical  indication  of  crack  growth  during  the  fracture  test  to  help  define  the 
compliance  pop-in  (Kj^,)  event.  An  uncalibrated  acoustical  indication  was  recorded  on  the 

Y  axis  simultaneously  with  load  on  the  X  axis  of  a  Houston  model  HR-97  X-Y  recorder.  The 
load  signal  was  obtained  by  attaching  an  auxiliary  slide  wire  to  the  load  indicator  of  the 
50,000-lb  tensile  machine.  The  slide  wire  was  one  leg  of  a  Wheatstone  bridge  measuring  cir¬ 
cuit.  This  gave  an  X  axis  record  of  the  bridge  unbalance  as  indicated  by  the  position  of  the 
wiper  arm  on  the  slide  wire.  The  bridge  output  was  calibrated  with  the  tensile  machine  load 
indication.  An  RCA  magnetic  cartridge  (1940  vintage)  with  a  0.050  inch  diameter  wire  pickup 
arm  inserted  in  the  notch  gave  an  indication  of  crack  activity.  The  acoustical  signal  was  am¬ 
plified  with  a  Ballantine  model  300  electronic  voltmeter,  filtered  with  a  20  KC  high  pass  filter 
to  eliminate  machine  noise,  and  rectified  for  the  Y  axis  input.  An  additional  monitor  was  ob¬ 
tained  by  using  a  high  impedance  Brush  crystal  head  set  in  parallel  to  the  X-Y  recorder.  A 
schematic  of  the  acoustical  monitoring  circuit  is  shown  in  Figure  4. 

Pictures  of  the  test  setup  and  a  close-up  of  the  compliance  gage  attached  to  the  specimen 
are  shown  in  Figures  5  and  6. 


PROCEDURE 

Wherever  applicable,  testproceduresusedwere  those  recommended  by  the  ASTM  Committee 
on  Fracture  Testing  of  High-Strength  Metallic  Materials  (References  3,  4  and  5). 

After  the  specimens  were  fatigue-cracked,  it  was  necessary  to  heat-tint  the  fatigue  crack 
surface  to  make  the  starter  crack  visible  for  optical  measurement.  Heat-tinting  was  done  at 
approximately  600°F  for  one  hour. 
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The  test  data  were  obtained  from  the  compliance  gage  recordings.  The  acoustical  indications, 
previously  described,  were  used  to  help  determine  the  pop-in  load  when  no  distinct  pop-in 
occurred  on  the  compliance  record.  In  some  cases,  the  deviation  from  linearity  of  the  com¬ 
pliance  slope  was  used  as  a  criteria  for  pop-in  load  determination.  The  lack  of  adequate 
elastic  constraint  in  the  sheet  materials  contributed  to  the  failure  to  obtain  distinct  pop-in 
load  values.  The  fracture-toughness  value  obtained  by  using  the  deviation  from  linearity  is 
commonly  designated  by  and  has  been  used  by  other  investigators  when  no  clearly  defined 

pop-in  has  occurred.  In  the  tabulated  data,  a  distinction  has  been  made  between  and 

For  some  of  the  materials,  the  acoustical  indications  helped  to  interpret  the  data;  however, 
the  use  of  acoustical  information  has  not  been  formally  endorsed  by  the  ASTM  Committee. 

A  thin  sheet  center  notch  specimen  has  a  predominate  plane  stress  state  rather  than  a  plane 
strain  state  required  for  conditions.  For  this  reason  the  center  notch  fracture-tou^iness 

specimen  is  not  considered  an  optimum  specimen  for  reproducible  type  measurements. 

Therefore,  the  K  data  presented  in  this  report  are  primarily  indicative  only  of  the  fracture 

initiation  resistance  properties  of  the  materials. 

Since  it  was  impossible  to  prevent  the  compliance  gage  knife  edges  from  moving  during  the 
latter  stages  of  crack  propagation,  optical  measurement  of  the  critical  crack  length  was  not 
conducted.  Instead,  the  sharp-notch  strength  has  been  used  to  indicate  the  fracture  resistance 
of  the  material.  The  sharp-notch  strength  of  a  material  is  defined  as  the  maximum  load  sus¬ 
tained  before  fracture  divided  by  the  nominal  cross-sectional  area,  i.e., 

P 

'^mox 

B  (W-20j,  ) 

where  P  is  the  maximum  load  sustained,  W  is  the  width,  B  is  the  thickness  and  2a  is  the 
max  ^ 

original  crack  length. 


COMPLIANCE  GAGE  CALIBRATION 

To  verify  the  compliance  calibration  curve  reported  by  Richard  W.  Boyle  (Reference  2), 
a  gage  calibration  test  for  center  notch  specimens  was  conducted  as  a  part  of  this  program. 

Nine  3  in.  by  12  in.  center-notch  specimens  of  AM  350  were  tested.  Simulated  cracks  were 
prepared  by  hand  sawing  slot  lengths  ranging  from  0  to  1.6  inches. 

In  this  method  of  compliance  gage  calibration  a  plot  of  Ev/crW  (ordinate)  vs.  vra/W  (abscissa) 
is  used  to  calculate  the  crack  length  (2a)  when  the  modulus  of  elasticity  (E),  specimen  dimen¬ 
sions,  and  the  deflections  at  all  loads  are  known.  Interested  readers  may  refer  to  Mr.  Boyle’s 
discussion  (Reference  2)  for  the  mechanics  of  using  the  calibration  curve. 

Figure  7  compares  our  calibration  curve  with  R.  W.  Boyle’s  curve.  The  excellent  corre¬ 
spondence  of  curve  shape  is  shown.  The  difference  of  ordinate  (Ev/crW)  values  is  accounted 
for  by  the  difference  in  Modulus  of  Elasticity  of  the  materials  used  in  the  tests,  i.e.,  steel 
versus  aluminum,  respectively. 

The  compliance  gage  calibration  curve  method  for  calculating  the  critical  crack  length  (2a) 
of  a  material  has  been  programmed  as  an  optional  part  of  the  computer  program  which  is  dis¬ 
cussed  later. 
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The  Irwin-Westergaard  analytical  basis  for  compliance  measurements  of  fracture  tough¬ 
ness  is  presented  in  Appendix  II. 


RESULTS  AND  DISCUSSION 

Table  2  shows  the  base-line  tensile  test  values  which  characterize  the  material  used  in  the 
fracture-tou^ess  tests.  The  yield  strength  value  for  PH  15-7  Mo  is  an  estimated  value  ob¬ 
tained  from  ejdsting  literature.  As  discussed  in  Reference  1,  an  error  of  10  percent  in  the 
yield  stress  will  result  in  a  2  percent  error  in  the  fracture-toughness  value  when  corrected 
for  the  plastic  zone  size. 

A  complete  tabulation  of  the  fracture  data  obtained  for  the  four  materials  can  be  found  in 

Tables  3,  4,  5  and  6.  The  initial  stress  intensity  value  of  K  has  been  tabulated  as  K  because 

nc 

the  majority  of  K  measurements  were  based  on  the  deviation  from  linearity  of  the  compliance 
gage  record.  Where  a  distinct  pop-in  occurred  it  is  noted  by  an  asterisk  in  the  data  tables. 

All  the  materials  tested  had  low  notch  strength  to  yield  strength  ratios.  The  following  tabu¬ 
lation  lists  the  fracture-tou^mess  values  obtained.  Gtenerally,  the  longitudinal  specimen  values 
were  higher  than  the  transverse  values.  A  longitudinal  specimen  has  its  major  axis  parallel 
to  the  rolling  direction  of  the  material. 


PH  15-7  Mo 
17-7  PH 
AM  350 
VJ  1000  (H-11) 


Longitudinal 

K 

nc 

Transverse 

Notch  Strength  Ratio 
Longitudinal  Transverse 

41.  3 

40.4 

.241 

.  264 

52.  6 

50.  3 

.798 

.  658 

60.  7 

54. 1 

.972 

.  891 

25.  3 

25.0 

.  142 

.  135 

Rockwell  C  hardness  readings  of  the  specimens  tested  showed  consistent  hardness  values 
over  the  entire  specimen.  The  Rc  hardness  ranges  were:  PH  15-7  Mo,  49.7  to  50.2;  17-7  PH, 
42.5  to  43.8;  AM  350,  44.5  to  45.5;  VJ  1000,  58.2  to  59.5. 

The  fracture  test  values  obtained  were  relatively  consistent  with  the  macroscopic  fracture 
surface  observations.  A  fracture  surface  with  granular  appearance  corresponds  to  cleavage, 
or  brittle  fracture,  while  a  fibrous  appearance  corresponds  to  a  ductile  failure. 


PH  15-7  Mo 
17-7  PH 
AM  350 
VJ  1000 


%  Shear  Appearance  of  Fracture 


Long. 

Trans. 

Longitudinal 

Transverse 

5 

5 

Rough 

and 

Granular  (G) 

85-90 

30-40 

Fibrous  (F) 

- 

Mixed  F  &  G 

95-100 

40-50 

Fibrous 

- 

Mixed  F  &  G 

5 

5 

Smooth 

and 

Granular 

The  difficulties  experienced  in  interpretation  of  some  of  the  compliance  curves  have  shown 
that  a  hi^er  magnification  than  the  550X  extensometer  used  in  this  program  is  desirable. 


Jones  and  Brown  (Reference  6)  have  recently  reviewed  the  use  of  acoustical  crack  growth 
detection  procedures  and  their  value  to  the  investigator.  Their  equipment  has  greater  sen¬ 
sitivity  and  sophistication  than  the  circuit  described  in  Figure  4;  however,  some  of  their  ob¬ 
servations  have  been  substantiated  by  this  program.  In  particular,  as  the  specimen  thickness 
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decreases,  the  signal-to-noise  ratio  decreases,  thus  permitting  extraneous  acoustical  indica¬ 
tions  prior  to  a  “pop-in”  or  a  deviation  from  linearity.  Extraneous  noise  is  a  problem  espe¬ 
cially  with  a  screw  driven  test  machine.  The  acoustical  pickup  serves  to  give  the  investigator 
an  additional  monitor  during  a  test  and  may  help  to  determine  the  “pop-in”  load  of  the 
material. 

The  only  comparative  data  values  found  in  the  literature  were  for  VJ  1000  in  Reference  4. 
A  mean  value  of  approximately  31  KSI  vHn.  for  found  therein  compares  favorably  with  the 

K,  value  of  25  KSI  yin.  and  the  K  values  from  30  to  33  KSI  obtained  from  the  com¬ 
ic  ^  c 

pliance  records  in  this  program. 


CONCLUSIONS 

1.  This  work  has  shown  that  of  the  alloys  tested,  VJ  1000  (H-11)  has  the  lowest  fracture 
toughness  while  AM  350  has  the  highest,  althou^  all  values  were  lower  than  expected  for  the 
particular  heat  treatments  tested. 

2.  The  longitudinal  and  transverse  fracture-toughness  values  were  very  similar,  although 
the  longitudinal  values  were  slightly  higher  for  the  materials  tested. 

3.  An  acoustical  monitor  is  of  value  in  determining  the  pop-in  load  where  a  definite  step 
pop-in  criterion  is  used. 

4.  Fracture-toughness  values  obtained  using  the  load  at  deviation  from  linearity  of  the  com¬ 
pliance  record  did  not  differ  appreciably  from  those  obtained  using  a  pop-in  load. 

5.  A  computer  program  greatly  simplifies  reduction  of  fracture-toughness  data,  especially 
when  the  compliance  method  is  used. 

6.  The  center-notched,  fracture-tou^iness  specimen  is  not  a  practical  specimen  for  deter¬ 
mination  of  plane  strain  fracture  toughness.  However,  it  is  highly  recommended  as  a  plane 
stress  specimen. 
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AM  350  (SCT  950)  VJ  1000  (2+2+2  at  lOOOP  F; 


Figure  1.  500X  Photomicrographs  of  Materials  Tested.  PH  15-7  Mo,  17-7  PH,  AM  350  were 
etched  with  a  solution  of  5g  FeCl3,  50  ml  HCl,  and  100  ml  H2O,  while  VJ  1000  was 
etched  with  4g  C6H2(N02)30H  and  100  ml  CH3OH2. 
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Figure  3.  Fatigue  Cracking  Setup  on  a  Six-Ton  Schenck  Fatigue  Machine 
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Z)  - 
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Figure  4.  Schematic  of  Acoustical  Indication  Circuit 


COMPLIANCE  GAGE 


AFML-TR-65-214 


11 


Figure  5.  Test  Equipment 
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Figure  6.  Close-up  of  Compliance  Gage  and  Acoustical  Pickup 
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Figure  7 .  Comparison  Compliance  Gage  Calibration  Curves 
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TABLE  2.  TENSILE  TEST  RESULTS 
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TABLE  3.  FRACTURE  TEST  RESULTS  FOR  PH  15-7  Mo 
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r*  2  Nominal  fracture  stress/yield  stress 


TABLE  3.  FRACTURE  TEST  RESULTS  FOR  PH  15-7  Mo  (cont. ) 
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TABLE  4,  FRACTURE  TEST  RESULTS  FOR  17-7  PH 
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TABLE  4,  FRACTURE  TEST  RESULTS  FOR  17-7  PH  (cont.  ) 
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TABLE  5.  FRACTURE  TEST  RESULTS  FOR  AM  350 
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TABLE  5,  FRACTURE  TEST  RESULTS  FOR  AM  350  (cont, ) 


AFML-TR-65-214 


20 


TABLE  6.  FRACTURE  TEST  RESULTS  FOR  VASCO  JET  1000 
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Indicates  occurrence  of  a  distinct  pop-in 


TABLE  6.  FRACTURE  TEST  RESULTS  FOR  VASCO  JET  1000  (cont.  ) 
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Indicates  occurrence  of  a  distinct  pop-in 
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APPENDIX  I 
COMPUTER  PROGRAM 


DISCUSSION 

The  computer  program,  described  below,  was  written  to  expedite  the  reduction  of  fracture 
toughness  ^ta.  Two  paths  of  calculation  are  available  as  distinguished  by  the  value  of  a  con¬ 
trol  variable  JJ.  When  JJ  is  set  equal  to  1  (one)  and  suitable  compliance  gage  data  are  sup¬ 
plied,  the  computer  calculates  the  critical  crack  length  of  the  material  before  calculating 

and  its  associated  parameters.  If  the  critical  crack  length  of  the  material  is  known,  the  com¬ 
pliance  calculation  will  be  bypassed  by  setting  JJ  equal  to  2.  The  data  input  format  for  JJ=1 
and  2  is  discussed  below. 

This  program  has  been  written  in  IBM  7094  Fortran  II  language.  To  provide  a  simplified 
means  of  data  input,  closed  subroutines  VDECOM  and  DECDCP  have  been  included.  When 
using  these  subroutines,  the  input  parameters  for  the  source  program  can  be  sequentially 
placed  on  data  cards  in  Columns  1  to  70.  However,  at  least  one  space  must  be  left  between 
each  entry  and  no  single  entry  may  overlap  to  the  following  data  card. 

The  critical  crack  length  calculation  path  (JJ=1)  has  been  programmed  using  the  Naval 
Research  Laboratory’s  procedure.  A  fifth  degree  polynomial  curve  fit  has  been  used  to  obtain 
the  following  equations: 


These  polynomial  coefficients  must  be  supplied  as  input  data  in  the  manner  described  below, 
regardless  of  the  computation  path  desired.  As  written,  the  compliance  calculation  requires 
four  sets  of  load-deflection  values  as  input.  The  first  three  sets  of  values  are  read  from  the 
linear  portion  of  the  compliance  curve.  These  values  are  averaged  to  obtain  a  representative 
value  of  Ev/crW  which  determines  the  magnitude  of  the  calibration  curve  shift.  The  fourth 
set  is  the  load-deflection  values  at  fracture. 

The  basic  fracture  toughness  calculations  are  made  using  Irwin’s  tangent  equation 

K  =  cr  W  font  - ) 

g  '  W  / 

The  corrected  K  values  incorporate  the  plastic  zone  correction  factor  where 

2  2  /  N 

a  =  a  +r„  and  K  Wtonf  -  +  - » — ) 

op  g  v  W  2W<7-  ^  / 

ys 

When  the  argument  of  the  tangent  fvinction  exceeds  the  mathematical  limitations  of  Irwin’s 
equation,  no  corrected  K  can  be  calculated,  thus  the  same  value  as  K  basic  is  printed  in  the 
data  output. 

When  all  the  input  parameters  necessary  to  make  fracture  toughness  calculations  are  known 
(JJ=2),  the  source  program  will  bypass  the  compliance  calculation  of  the  critical  crack  length. 

An  example  of  the  data  output  from  this  computer  program  is  shown  in  Table  7. 
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DATA  INPUT 

All  input  data,  except  JJ,  must  be  decimal  values.  The  first  two  data  cards  for  both  JJ=1 
and  JJ=2  contain  the  polynominal  curve  fit  constants  shown  below.  The  minimum  one  space 
format  using  columns  1  to  70  apply  to  all  the  data  card  input  discussions.  The  constant  data 
cards  numbers  1  and  2  are  supplied  only  once  for  all  program  runs. 

Card  No.  1 

.3193?  .£2331  .15115  0.  0.  .38496 

•00000000000000000000000  000  000000000000000000000000000000000000000000000000000 

I  2  3  4  5  i  7  B  9  10  11  12  13  14  15  16  17  11  ts  30  21  22  23  24  2S  28  27  2S  23  30  3t  32  33  34  35  36  37  38  39  40  41  42  43  44  45  46  47  48  4B  50  SI  U  53  54  55  56  97  51  39  60  fil  62  63  M  65  SB  87  BB  SB  70  7t  72  73  74  75  76  77  71  7B  V 

II  1111111111  111  1  11111111111111111111111111111111111111111111111111111111111 

222222222  222222222222222222222222222222222222222222222222222222222222222222222 
333333  333  333  33333333  333  33  3333333333333333333333333333333333333333333333 

44444444444444444444444444444444444  44444444444444444444444444444444444444444444 
SS5S55555535555505  55  5555555555555555555555555555555555555555555555555555555555 
6t666S66S6B666S666866666666668666E666  888888888888888888888888888888888888888888 
7  7  7  7.7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  /  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7 
881  888  88  8888  88888888  888  88  8  888888883888888888888888888888888888888788888 
999  9  999999999999999999999999999999  9999999999999999999999999999999999999999999 

I  2  3  4  5  E  7  8  S  10  II  12  13  14  15  IB  17  H  If  a  21  22  23  24  25  2S  27  28  23  30  31  32  33  34  35  36  37  38  39  0  41  42  43  44  45  46  47  40  40  50  51  52  53  54  55  56  57  5!  59  60  El  62  E3  64  65  66  67  68  63  70  71  72  73  74  75  76  77  78  79  80 

'T  A  83  O  O  ^  Q  I 

Hi.ni  SO  0  1 


Card  No.  2 

1.E3134  -.£1£053  -.0774344  0.  .0019433 


00080000000000  000000  00000000  0000  0000000000000000000000000000000000000000000 
1  2  3  4  5  I  7  I  9IOIII2l3l4l$l6t7l8l9  20  21  22  23  34S  26  27  2l29  30  31  32  33  34  3S36  37  3l3S4a4l«7  43  44  4S  46  47  48  49SOSIS?S3S4S5  56  $7S8S9  60  $l  62  63  64  65  G6  67  S8  69  70  7  l  72  73  74  75  78  77  7|76BO 
111  1111111  111111111111111111111111  111111111111111111111111111111111111111111 

22  22222222  2  22222222222  2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 

3  3  333333  33333  333  3333333333  33  333333  33333333333333333333333333333333333333 

444444  44444444444444444  4  44444444444  4444444444444444444444444444444444444444 

555555555555555  5555555555555555555555555555555555555555555555555555555555555555 

8BSS5BE88BB8E86858BE858EEE88B8B58BE655B88BEBE8E58B85BB8888BEEEEE8B85B8BE5E8EEB8E 


77777777777777777777/7 
8  888  8888  888888888  8888888888  88  88888  883888888888888888888888888688888388888 
99999999999999999999999999999999999999  99999999999999999999999999999999999999999 

1  2  3  4  J  soe^'^  '*  '5  22  n  24  25  H  27  a  26  30  31  32  33  34  35  36  37  31  141  42  43  «  45  46  47  «  49  50  51  53  53  54  55  56  5?  M  59  60  61  62  63  61  65  66  61  66  65  70  II  72  13  74  15  76  77  71  79  to 
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JJ=1:  Actual  test  data  is  read  in  on  subsequent  cards,  beginning  with  card  number  3.  Each 
set  of  test  data  is  begun  in  column  one  and  may  be  presented  on  several  successive  cards,  as 
necessary.  A  value  must  be  given  for  every  parameter  listed  on  the  card  below.  If  applicable, 
0.0  may  be  inserted  as  the  value  for  any  of  these  items. 

Sample  Data  Card  for  JJ=1 


-e- 

< 

I 


$  OD 

I  I 


000000000000000  oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

1  2  3  4  5  6  7  8  9  10  11  12  »3  14  15  16  17  18  19  20  2l  22  23  24  25  26  27  28  29  30  31  32  33  34  35  36  37  38  39  40  4)  42  43  44  45  46  47  48  49  50  51  52  53  54  55  56  57  58  59  60  61  62  S3  64  65  66  67  68  69  70  71  17  73  74  75  76  77  78  79  80 

' ' '  II  1 1 1 1 1 1 1  1 1  I  I  I  I  1 1 1 1 11 1 1 11 1 n 1 1 11  n  1 1  1 1 1 1 1 1 1 1 1 1  I  111111111111111111 
22222222222222222222222222222222  222222222222222222222222  22  2222222222222222222 
33  333333  3333  3333333  33333  3333  3  333  3  333  33  333  333  333  33  33333333333333 

444444444444444444444444444444444444444444444444444  4444444444444444444444444444 
55555555555  555555555555555555  555555555555555555555555  555555555555555555555555 
066  66666666  666666666666666666666666666666666  666666666666666666666666666666666 
77777777777777777  77777777777777777777  77777  777777/7777777777777777777777777777 
88  888888  8888  8  88888  888  8  88  8  88  8888  8  888  88  88  888  888  8  88888888888888 
9  9  9  9  9  9  '  9  9  9  '  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9 

I  2  3  4  5  6  7  8  9  10  11  12  13  14  15  16  17  18  19  20  21  22  23  24  25  26  27  28  29  30  31  32  33  34  35  36  37  39  39  40  41  42  43  44  45  4f.  47  43  49  SO  51  52  53  54  55  56  57  58  59  60  61  62  63  64  65  66  67  68  69  70  71  72  73  74  75  76  77  71  79  10 
5081  esc 


NOTE:  This  data  card  yields  the  first  two  lines  of  results  in  Table  7. 

JJ=2:  For  calculations  not  requiring  the  compliance  calculation  of  the  critical  crack  length, 
the  following  sample  data  card  is  appropriate.  The  one  space  between  data  items  in  Columns  1 
to  70  is  required. 
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Sample  Data  Card  for  JJ=2 


> 

u. 

1 

< 

1 

1 

CD  q! 

1  1 

to 

< 

Q. 

t- 

.  936 

1 

3 .014 

1  1 

, 052  4.45 

1 

225. 

1.91 

1 

7.  75 

OOOOGOOOO  0000  00000000000000000000000000000000000000000000000000000000000000000 

1  2  3  4  5  6  7  9  9  >0  II  12  13  14  15  18  17  18  19  M  21  22  23  24  25  26  27  28  29  30  31  32  33  34  35  36  37  38  39  40  41  42  43  44  45  46  4I  44  49  50  SI  U  53  54  55  56  57  51  59  60  81  62  63  M  65  S6  67  M  69  70  71  72  73  74  75  76  77  7|  7|  IQ 

1 1 ' ' '  I ' '  >  >  <  1  n  1 11 1 1 1 1 1 1 1 11 1  11  111111111111111111111111111111111111111111111111 
222222222222222  222222  2222222222222222222222222222222222222222222222222222222 
33  3333  3333  33333  333333  33  3333  333333333333333333333333333333333333333333333 
44444444444  444444  4  44444444444444444444444444444444444444444444444444444444444 
5 1  5  5  5  5  5  5  0  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5 
60066  66666666666666666666666666666666686666666668666686666666666666666666666666 
mmmmnmmnmiimm  i  m  i  m  m  i  m  mm  m  m  m  m  n  n  n  m  m 
88  a  838  8888  88888  888888  88  8888  888888888888888888888888888888888888888488888 
399  99939999999999999999999999  9999999999999999999999999999999999999999999999999 

1  2  3  i  5  6  7  8  1  to  11  12  13  14  55  IS  17  13  1?  ?0  21  23  23  J-i  25  2*  ?7  20  29  30  31  32  33  3<  35  36  37  38  3?  0  41  42  43  44  45  46  47  48  49  50  51  52  53  54  55  56  57  58  59  50  $I  62  63  64  65  S$  67  68  69  70  7l  72  73  74  75  76  77  78  79  M 
soei 


The  following  is  an  abbreviated  glossary  for  the  symbols  shown  on  the  sample  data  cards. 
No  attempt  has  been  made  to  define  all  the  symbols  used  in  the  source  program  and  the  sub¬ 
routines  VDECOM  and  DECDCP. 


Fortran  Symbol 
A(p 
W 
B 

PIC 

SY 

FAC 

E 

P(1),P(2),P(3) 

DEF(l),  DEF(2),  DEF(3) 

P(4)  or  PF 
DEF  (4) 

A 


LIST  OF  SYMBOLS  (INPUT  DATA) 

Analytical  Symbol  Definition 

2a.Q  Total  initial  crack  length 

w  Width 

Thickness 

"Pop-in"  load  (Kips) 

Yield  Stress  (I^I) 

Total  magnification  of  the 
compliance  gage* 

E  Modulus  of  Elasticity  (psi) 

Loads  (Kips)  on  linear 

portion  of  compliance  curve 
Deflections  (inches  of  chart) 

corresponding  to  P(l),  P(2),  P(3) 
Load  at  Fracture  (Kips) 

Deflection  at  fracture  (inches 
of  chart) 

2a  Total  critical  crack  length 


*  A  gage  with  a  250X  microformer  and  a  gage  lever  ratio  of  2  would  have  FAC  =  500. 
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DATA  OUTPUT 


The  computer  program  supplies  a  data  output  as  shown  in  Table  7.  The  first  line  of  a  single 

test  set  contains  K_  data  while  the  second  line  contains  K  data.  The  “nominal  ratio”  as  used 
Tc  c 

here  is  equivalent  to  the  “notch  strength  ratio.” 


SOURCE  PROGRAM  AND  SUBROUTINES 


The  complete  Fortran  11  program  which  consists  of  the  source  program  and  two  closed  sub¬ 
routines,  VDECOM  and  DECDCP,  is  presented  in  Table  8. 
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TABLE  8.  SOURCE  COMPUTER  PROGRAM 

CCALC  CENTFR  NOTCH  FRACTURF  TOUGHNESS  AND  COMPLIANCE  GAGE  CALCULATIONS 
DIMENSION  AFIT(7)»  RFIT(7)»  P(4)»  V(S)»  DEE(4),  Y(3)»  INT{10)» 

I  DFC(IO) 

PI  =  3.14159 
JTAPP  =  2 
NTAPF  =  3 

READ  IN  COEFFICIENTS  FOR  COMPLIANCE  CURVES 
Y  =  F(X  ) 

READ  INPUT  TAPE  JTAPF,  1000.  (  AFIT(I).  I  =  1.  7) 

C  X  =  F( 1  ./Y  ) 

READ  INPUT  TAPE  JTAPF,  1000,  (  BFIT(I)  ♦  I  =  1*  7) 

WRITF  OUTPUT  TAPE  NTAPF, 4000 
5  N  =  1 

KPASS  =  1 

CALL  VDECOM(  M,  INT,  D^^C,  KPASS) 

JJ  =  INT (  1  ) 

N  =  7 
KPASS  =  3 

call  VDFCOMI  M,  ImT,  DFC.  KPASS) 

AO  =  PFr( 1 ) 

AO  =  AO 
W  =  nFC(2) 

B  =  OFC { 3 ) 

Die  =  DFC (4 ) 

SY  =  DFC ( 5  ) 

FAC  =  DFC(6) 

F  =  DFC(7) 

A  =  DEC(6) 

PF  =  Pf"C(7) 

MM  =  0 

GO  TO  (  1  5 , 35  )  ,  JJ 
1  5  M  =  5 

KPASS  =  3 

call  VDEC0M(  N,  INT,  DEC,  KPASS) 

PO  10  I  =  1»  4 
K  =  2  «■  I 

P ( I )  =  1000. *DEC ( K-1 ) 

1  0  DEF ( I )  =  DEC ( K ) 

C  COMPUTE  AVERAGE"  VALUES  nc  y  AT  7 

YSUM=  0. 

DO  20  I  =  1  .  S 

V( T )  =  0.5  #  OFF ( I )  /  fat 

Y(T)=  E^^V(I)*R/P(T) 

7  0  ycum^tYSUM  +  Y(I) 

yavg=ysum/s. 

YFRAC=R*B»DRF ( 4 ) /{2.*FAr*P{4) ) 

7=pT*A0/(  2.  *'>"') 
r  FVALUATE  COMPLIANC*^  POLYNOMIAL  AT  ? 

YF  IT  =  AF  IT ( 1 ) 
no  25  I  =  1,  5 

25  YFIT  =  YFIT  +  AEIT(I+1)  *  Z**I 
C  COMPUTE  DIFFFRFNCE  RETWFFN  Y  VALUI^S 

YOIFR  =  YAVG  -  YFIT 

C  CALCULATE  AT  FRACTURE  FROM  ADJUSTED  VALUE  OF  Y 

YNFW=YFRAC  -  ydtfr 
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TABLE  8.  (Continued) 

7MFW  =  RF  I  T (  1  ) 
no  RO  I  =  T  »  ^ 

•:<n  7NFW  =  7NFW  +  BFIT  (I  +  l)  /(yn^W**!) 

C  COIYDUTF  NEW  PRACTURF  CRACK  LFMGTH 

A  =  ?,  *  W  *  7NPW  /  PI 
DF  =  P( 4 ) /I  000. 

C  C0Y1PUTF  THE  STRESS  AND  pR  AC  T(  IR E-TOUGHNF  SS 

RE  AKIC  =  0.0 

EG  =  PIC/(W*R) 

FNP  =  PIC/(R^f  (W-AO)  \ 

FM1  =  PMP 

AN?  =  DIC/(R^f  (W-Al  )  ) 

FRN  =  GNP/GV 
PA  =  PI*A0/{2.*W) 

RKP  =  GI  NF  (  PA  ) /rO.PF  (  D  A  )  ^fW*PG**2 
AOW  =  AO/W 
SR  =  SNF/SY 

RP  =  RKS/ ( 2.*pI*Sy**7 ) 

RKl  =  SQRTFIRKS) 

C  evaluate  the  EOUAIION  LIYilTATION 

FOYP  =  (SG/SY1**2 
QM  =  SORTF ( 2 . /GOYS  -  1.) 

PAY  =  ATANF(OV)  -  Ovi*gOYc/2. 

I F ( PA V-DA  )  RR  ,40  »  40 
RQ  AKIC  =  RKl 
GO  TO  45 

C  COY'PUTF  THP  PLASTIC  70MF  FR ACTURF-TOUGHNpSS 

40  PAP  =  PA  +  RKC/ ( 2.0*W*SY*»2 ) 

PKP  =  S INF { PAP ) /COSF ( PAP ) *W*PG**2 
RK  =  SORTF(RKP) 

AKIC  =  SQRTF(PKS) 

RKS  =  PKS 

IF ( AK I C-RK-. 005 ) 44,40,40 

44  RP  =  PKS/ (  2  ,^fP  P^FY**2  ) 

45  I F ( MM ) 50 . 50 .60 

50  WRITE  OUTPUT  TAPE  N T AP E , 2000 ♦ AO » W , R , P I C , SG * SN 1 . SNF , SY , SR , AOW » RP » 

1  Ri^  1  ,AK  IP 
A1  =  AO 
AO  =  A 
PIC  =  PF 
NN  =  1 

C  RFPFAT  FOR  PlANF-STRESS  TOUGHNESS 

I F ( AO ) 5 ,5  ,R5 

60  WRITE  OUTPUT  TAPE  NT APE » ROOD . A  1  ♦ AO , W , R ♦ P I C , SG . SN2 . SNF , SY , SRN ♦ SR  , 

1  RP»RK1  ,AKIC 
r-0  TO  c, 

1000  FORMAT!  7F10.5) 

2000  FORMAT ( /,F7.R,F15,R,F9.4,F8.2»F7.2,F9.2,F8.2.F7.2,F14,2,F8,2.F9.4. 
1  2F9,2) 

ROnO  format (  F7.R,F8.R»F7.R,F9.4,F8.2,F7.2.F9.2.F8.2.RF7.2,F17.4. 

1  2F9.2) 

4000  FORMAT ( 1  HI  ,  1 16H  CRACK  LENGTH  WIPTH  THICKNESS  LOAD  GROSS  NOM 

IINAL  NET  YIELD  NOMINAL  NET  LENGTH/  PLASTIC  TOUGHNESS  K, 
2/,120H  INITIAL  FINAL  STRESS  STRESS  STR 

RESS  STRESS  RATIO  RATIO  WIDTH  70NE  BASIC  CORRECTED) 

END 
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TABLE  8.  (Continued) 


SUBROUTINE  VDECOMIN 


KARRAY.P  ,KPASS)  12/lfl/t3 


C 

c 


c 

B 

B 

B 

B 

B 

B 

C 


C 


C 


C 


C 


SUBROUTINE  VDECOMN  ,  KARRAY.P  .KPASS) 

dimension  k1nPUT(72),  KARRAY(IO),  PARRAY(10),AlNPLr(72), 
I  KCUTPT{72).  ACUTPT172) 

equivalence  ( plus .  IPLUS ) . { am INUS. I M 1 NLS )  ,  (DECPT ,  IDECPT  ) 
EUUivALfcNut  ip),(  blank.  IBLANK). 

1  - 


I  /\  T  MDl  T  -  K  I  MPIj  T  ) 


i  ADI  TOT 


KPASS  =  KPASS 
J2  =  INPUT  TAPE 
J2  =  2 

J4  =  SCRATCH  TAPE 
JA  =  A 

GC  TO  (  2,  2.  2,  3),  KPASS 

2  REWIND  JA 

3  NLMCCP  =  N 
NEX  =1 

N1  =  I 
L  =  I 

GC  TO  (I.  1C.  5C.  50).  KPASS 
1  CONTINUE 

SET  UP  characters  FOR  LATER  TEST 
PLUS  =  206060606060 
AMINUS  =  A06060606060 
OECPT  =  336060606060 
COMMA  =  736060606060 
E  =  2560606C6C60 

blank  =  6C6060606060 

READ  ALPHANUMERIC  CHARACTERS 
5  READ  INPUT  TAPEJ2.100C.  ( A  I NPUT ( J  )  ,  J  =  1 .  72  ) 


1  =  1 

GC  TO  (1C. 1C.  50,  50),  KPASS 
DECOMPCSITION  OF  INTEGERS 

10  DO  21  N  =  Nl,  NUMDCP 
N1  =  N 

SEARCH  FCR  START  OF  NUMBER 
ICl  IF  (KINPUT(I)  -  IBLANK)  102.  Ilf  102 
102  IF  (KINPUT(I)  -  ICOMMA)  12  .  11.  12 

11  1=  I+l 

1F(  I  -72)101,101,  5 
SELECT  INTEGERS 

12  LI  =  L 
M  =  0 

DC  20  J=l,6 

KOUTPT ( L )  =  K  INPUT { I  ) 

IF(  KINPUT(I)  -  IMINUS)  122,  120.  122 

122  IF(  KINPUTII)  -  IPLUS)  123,  120,  123 

120  M  =  1 

123  L  =  L  +  1 
1=  I  +  l 

IF(KINPUT(I)  -IBLANK)121  ,13  ,121 

121  IF(KINPUT(I)  -ICOMMA)  20  ,13  ,  20 

RIGHT  ADJUST  IN  FIELD 

13  IF(  J-6)  1A,21,14 

lA  KDO  =  J  -  M 

DC  15  K  =  1,  KDC 
L2=L1+6-K 
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TABLE  8.  (Continued) 


SUBROUTINE  VDECOMIN  ,  KARRAY.P  ,KPASS) 

L3  =  LI  +  J  -  K 

15  KCUTPTI  L2  )=  KCUTPK  L3  ) 

L  =  LI  ♦  6 

KDO  =  6-J  +  P 

KGO  =  1  +  M 
DC  16  K  =  KGOf  KDC 
L4  =  LI  +  K  -1 

16  KCUTPK  L4)  =  0 

GC  TO  21 

20  CCNTINUE 

21  CCNTINUE 

lEND  =  6»  NLMCCP 

C  WRITE  ALPHANUMERIC  CHARACTERS 

WRITE  OUTPUT  TAPEJA.IOOO,  (AOUTPT(J),  J=lf  lEND) 
REWINCJA 

C  READ  INTEGER  LIST 

READ  INPUT  TAPEJ4,1001,  (KARRAY(J),  J=1,NUMDCP) 
REWINDJ4 
AO  RETURN 

50  CALL  CECDCP  (  IPLUSt  IMINUSt  ICECPT,  ICOMMA, 

1  KINPUTt  NUMDCPt  PARRAYt  It  Lt  NEXt  Nit  KPASS) 

NEX  =  NEX 
GO  T0(30,5),  NEX 

30  DC  31  J  =  It  NUMDCP 

31  P(J)  =  PARRAY(J) 

GC  TO  AC 

1000  FCRMAT(72A1) 
lOCl  FCRMAT(11I6) 

ENO( I,l,0,0,0,ltltlf0,lf0,0,0,0,0) 


12/18/63 


IE,  lELANK, 
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TABLE  8.  (Continued) 


SUBROUTINE 

DECDCPI 

IPLUS, 

IMINUS, 

IDECPT 

,  ICOMMA, 

IE, 

IBLANK, 

12/18/63 

SUBROUTINE 

CECDCP( 

IPLUS, 

IMINUS, 

IDECPT 

,  ICCMMA, 

IE, 

IBLANK, 

KINPUT,  NUMDCP,  PARRAY, 

I,  L,  NEX 

,  Nl, 

KPASS ) 

EQUIVALENCE  (AOUTPT,  KOUTPT) 

DIMENSION  AOLTPT(72),  KOUTPT(72),  PARRAY(IO)  ,  KINPUT{72) 
NEX  =  NEX 

C  J2  =  INPUT  TAPE 

J2  =  2 

C  J4  =  SCRATCH  TAPE 

=  4 
M  =  L 

GC  TO  (50, 51), NEX 

DECOMPOSITION  OF  DECIMAL  AND  EXPCNENTIAL  NUMBERS 
LIMIT  DECOMPOSITION  TO  6  NUMBERS 

50  IFINUMDCP  -6)503,503,502 

502  IDEC  =  6 
lEND  =  12«ICEC 
GC  TO  50A 

503  IDEC  =  NUMDCP 
lEND  =  12»  NUMDCP 

5C4  IF{ I  -72)  51,51,505 
505  NEX  =  2 

GO  rc  300 

51  NEX  =  1 
DO  IOC  N=  M,  IDEC 
Nl  =  N 

C  SEARCH  FCR  START  OF  NUMBER 

510  1F(KINPUT(  I  )-  IBLANK)  52,  53,  52 

52  IF (K  INPUT (  I  )-  ICOMMA)  54,  53,  54 

53  1=  I  +  1 

IF(  I  -  72)  510,  510,  505 

54  Ml  =  N 

C  STORF  NUMBERS  UP  TO  DECIMAL  POINT 

541  1F(KINPUT(I)  -  IDECPT)  55,  65,  55 

55  KOUTPT (M )  =  K  INPUT ( I  ) 


1=  I  +1 
M=  M  +1 
GC  TO  541 


TEST  FOR  END 

OF  NUMBER 

OR 

EXPCNENTIAL 

60 

1  F ( K  INPUT (  I  )  - 

IE 

) 

61, 

70,61 

61 

IF(K  INPUT!  I  )  - 

IPLUS 

) 

62, 

70,62 

62 

IF(K1NPUT(  I  )  - 

IMINUS 

) 

63, 

70,63 

63 

IF(KINPUT(  I  )  - 

ICCMMA 

) 

64, 

80,64 

64 

IF(KINPUT(  I  )  - 

IBLANK 

) 

65, 

80,65 

STORE  DECIMAL 

POINT 

AND 

NUMBERS 

65 

KCUTPT (M  )  =  K  INPUT ( I  ) 

1=  1*1 
M=  M+1 
GC  TO  60 

C  COMPLETE  EXPCNENTIAL  FIELD  THROUGH  8  LCCATICNS 

70  LDO  =  Ml  ♦  7 

DC  71  Jl=  M,  LDO 

71  KCUTPT(Jl)  =  0 
M  =  Ml  ♦  8 

C  STORE  E  IN  LOCATION  9 

KCUTPKM)  =  IE 
IF(KINPUT( I )-  IE)  73,72,73 
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TABLE  8.  (Continued) 


SUBROUTINE  CECDCPI  IPLUS,  IWINUS,  ICECPT,  ICONM«,  IE,  IBLANK, 


72 

1=  I  +  l 

TEST  FCR  SIGN  OF  EXPONENT 

73 

IFIKINPUTI I )- 

IMINUS)  7A, 

76, 

7A 

7A 

IFIKINPUTI I )- 

IPLUS  )  75, 

76 

,75 

75 

KCUTPT(M+1)  = 

IPLUS 

GC  TO  77 

STORE  SIGN 

76 

KCUTPT(M+1)  = 

KINPUTI  I  ) 

1  =  1  +  1 

TEST  FCR  END 

OF  EXPONENT 

77 

IFIKINPUT (  I  +  l )  -  I  BLANK) 

78, 

79,  78 

78 

IFIKINPUTI I+l )  -  ICOMMA) 

791 

,79  ,791 

79 

KCUTPTIM+2)  = 

0 

KCUTPTIM+3)  = 

KINPUT I  I ) 

GC  TO  792 

791 

KCUTPTIM+2)  = 

KINPUTI I ) 

KCUTPTIM+3)  = 

KINPUTI  I  +  l) 

792 

1=  1  +  2 

M  =  Ml  +  12 
GC  TO  ICC 

C  COMPLETE  DECIMAL  FIELD 

80  LDO  =  Ml  +  11 

DO  81  J1  =  M  ,  LDO 

81  KCUTPTI J1  )  =  0 
M  =  Ml  +  12 

ICO  CONTINUE 

C  WRITE  ALPHANUMERIC  CHARACTERS 

WRITE  OUTPUT  TAPE J A ,  1 000 , I ACUT PT { J ) , J=  1 ,  lEND) 
IFINUMDCP  -  ICEC  )  201,  201,  200 

200  M  =  1 

N1  =  N1  +  1 

lEND  =  12*(NUMDCP  -  6) 

ICEC  =  NUMDCP 
GC  TO  50A 

201  GC  TUI  202,  2C2,  202,  300),  KPASS 
2C2  REWINDJA 

C  READ  DECIMAL  AND  EXPONENTIAL  LIST 

READ  INPUT  TAPE  J A , 1 0 10 , ( P ARRA Y I J ) , J =  1  ,  NUMDCP) 
REWINDJA 
300  L  =  M 
RETURN 

1000  FCRMAT(72A1) 

1010  FCRMAT(6F12.5) 

END!  1,1, 0,0, 0,1, 1,1, 0,1, 0,0, 0,0,0) 


12/18/63 
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APPENDIX  II 

COMPLIANCE  METHOD  ANALYSIS 


From  the  theory  of  elasticity  as  presented  by  Timoshenko  and  Goodier  (Reference  7) ,  there 

a4 

is  a  stress  fimctionc^  which  satisfies  the  biharmonic  equation  A  c^=0  and  also  satisfies  bound¬ 
ary  conditions  such  that  the  external  forces  may  be  considered  as  an  extension  of  the  internal 
stress  distribution.  If  this  is  true  then: 


d  ^ 
dy* 


ond  T 


xy 


dx  dy 


where  is  a  function  of  x  and  y  and  is  called  the  Airy  stress  function. 


Westergaard  (Reference  8)  has  provided  a  two-dimensional  stress  analysis  of  a  very  large 
flat  plate  with  tension  applied  in  the  y  direction  and  a  system  of  cracks  along  the  x  axis  each 
of  length  2a  and  center  atx  =  o,tL,±.2L —  Using  Westergaard’s  notation,  if  the  function  is 
related  to  a  complex  function,  we  shall  call  it  Z  where 


Z  =  Z(z)=Z(x+ly)ReZ-t-iImZ  (I) 

and  the  pertinent  values  are  given  by  the  analytic  function 


In  Westergaard’s  notation. 


d  z 


-  dz 

z  =  - 

d  z 


Applying  the  Cauchy  Riemann  conditions. 


dRe  Z 

a* 


aim  Z 

ay 


R8  Z  ,  and 


dimZ 

a* 


aReZ 

ay 


Im  Z' 
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If  the  Airy  function  is  defined  as 

<^  =  Re  T -I-  y  im  z 


and  it  satisfies  the  biharmonic  equation 

4  4 

d  4>  zd  4> 

<3 


then 


O  Re  Z 

ax2 


0 


y 


d  Im  z 
ax2 


Since 


and 


Therefore, 


a  Re  Z 

ax^ 


=  Re  Z 


a  Im  Z 

~d7^ 


=  y  Im  Z 


o-y  =  Re  Z  +  y  Im  z' 


By  similar  steps, 


=ReZ-yrmZ 


T 


xy 


-y  Re  z' 


(3) 


(4) 

(5) 


The  displacement,  v,  in  the  y  direction  for  a  plane  strain  situation  is  given  by  the  equation. 

Ev  =  2  (l-vZ  )rm  Z-(l4v  )y  ReZ 
where  v  is  Poisson’s  ratio  and  E  is  the  modulus  of  elasticity. 

It  has  been  shown  by  Irwin  (Reference  9)  that  the  above  stress  analysis  is  a  good  approxi¬ 
mation  to  the  stress  state  in  a  centrally  cracked  sheet  specimen  whose  vertical  and  horizontal 
axes  of  symmetry  are  taken  as  they  and  x  axes  respectively.  The  specimen,  with  width  w  =  L, 
is  regarded  as  one  unit  of  the  crack  system,  and  therefore  the  specimen  edges  are  at  x  =  ±L/2 
and  the  crack  extends  from  x  =  -a  to  x  =  +a.  Boundary  conditions  require  the  stresses  a  and 

^  zero  along  the  borders  of  the  crack  and  the  stresses  c  and  r  to  be  zero  along 

X  xy  ° 

the  side  boundries.  All  these  conditions  are  fulfilled  except  the  condition  that  cr  =0.  Irwin 
remedied  this  by  rewriting  c  as  ^ 


<r  =  Re  Z  -y  I  m  Z  '  “ 

X  ox 


(6) 
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where  a 


o 


is  a  constant  stress  adjusted  so  that 


X 


0  at  X  =  ±  — 
2 


o 


The  constant  a  does  not  make  a  zero  at  all  points  on  the  side  boundries,  but  it  does  reduce 

O  X 


cr^  to  a  very  small  value.  It  is  assumed  here  that  this  deviation  from  an  exact  solution  results 

in  errors  much  smaller  than  those  resulting  from  the  departure  from  linear  elasticity  theory 
for  finite  strains. 

Irwin  modified  the  equation  for  the  displacement  in  the  y  direction  to  make  it  applicable  to 
a  generalized  plane  stress  condition.  This  equation  is 


(7) 


Ev=2rfnZ  —  (l-J-vjyRe  Z 


It  is  now  possible  to  derive  an  expression  for  compliance  at  the  specimen  edges.  The  edges 
of  the  specimen  are  at  x  =  ±L/2.  Therefore,  z  =  ±L/2  +  iy  where  y  is  1/2  of  the  specimen 
gage  length 

Integration  of  Z  =  f  Zdz  utilizing  Equation  (2)  gives 


(8) 


where  C  is  a  complex  constant  of  integration  to  be  evaluated  for  the  known  conditions  at  a  =  o. 
Utilizing  the  identities 


cos 


and 


and  substituting  in  Equation  8,  it  is  found  that 


+  Re  C  +  i  Xm  C 


(9) 


cos 


L 
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Therefore, 


im  z 


Also,  when  x  =  ±L/2, 


L  tr 

TT 


sinh 


sinh 


cos 


^■y 


L 


Im  C 


(10) 


(II) 


and  it  is  real.  Therefore,  Re  Z  =  Z  at  x  =  ±L/2.  The  expression  for  tlie  displacement  v  is 
found  by  substituting  Equations  10  and  11  into  Equation  7.  In  accordance  with  R.  W.  Boyle 
(Reference  2) ,  the  compliance  v  /c ,  is  expressed  as  a  dimensionless  parameter  by  dividing 
by  the  specimen  width  w  and  multiplying  by  Young’s  modulus.  Now,  by  setting  L  =  w,  the 
equation  for  compliance  becomes 


Ev 

sinh 

^y  ^ 

W 

\  ,  ZImZ 

y 

<r  W 

1 

1  COS 

TTQ 

w  y 

j  aw 

W 

( I  +v  ) 


(in^ 


I  - 


cosh 


2 


■ra 

w 

LL 

w 


The  term  Im  C  is  evaluated  by  knowing  that  when  a  =  o  the  specimen  is  simply  an  uimotched 
specimen,  and  E  v  /cr  W  is  equal  to  y/w.  The  final  equation 


Ev 

(T  W 


—  sinh 

TT 


(  n-  V  ) 


<  y 

w 


(12) 


To  review  the  pertinent  symbols  in  Equation  12.  E  is  Young’s  modulus,  xis  one-half  the  speci¬ 
men  extension,  a  is  the  gross  stress  calculated  from  the  dimensions  of  the  uncracked  speci¬ 
men,  _^is  the  specimen  width,  ^  is  one-half  the  crack  length,  and^is  one-half  the  specimen 
gage  length. 

It  can  be  seen  from  Equation  12  that  a  plot  of  the  dimensionless  parameters  Ev  /aW  vs 

will  yield  a  curve  whose  shape  is  constant  regardless  of  the  material  tested.  This  is  true 

because  the  factor  y/W  is  constant  for  identical  specimen  geometries.  The  curve  will  shift 
up  or  down  the  Ev  /crW  axis  for  different  materials  because  of  different  E  values  and  the 
corresponding  change  in  the  extension  factor  vy/W. 
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